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Abstract 


In this paper the elasticity problem for a thick plate with a circular 
hole is formulated in a systematic fashion by using the z-component of 
the Galerkin vector and that of Muki's harmonic vector function. The 
problem was originally solved byAlblas [1]. The reasons for reconsidering 
it are to develop a technique which may be used in solving the elasticity 
problem for a multilayered plate and to verify and extend the results given 
by Alblas. As in [1] the problem is reduced to an infinite system of 
algebraic equations which is solved by the method of reduction. Various 
stress components are tabulated as functions of a/h, z/h, r/a, and v, 
a and 2h being the radius of the hole and the plate thickness and v the 
Poisson's ratio. Among the addi tonal results of particular interest one 
may mention the significant effect of the Poisson's ratio on the behavior 
and the magnitude of the stresses. 


1. Introduction 

In this paper the three-dimensional elasticity problem of a thick 
plate containing a circular hole and subjected to unidirectional loading 
away from the hole region is reconsidered. The problem was originally 
solved byAlblas for stretching as well as bending of the plate [1]. Alblas 
obtained the elasticity solution of the problem by formulating it in terms 
of three displacement potentials and by reducing it to an infinite system 
of linear algebraic equations. Basically this system of equations arises 
from the expansion of the boundary conditions on the hole surface into 
series of complex eigenfunctions. The infinite system was theh”¥6l"ved 
by the method of reduction and the results were compared with varous existing 
approximate solutions. 

The main reason for reconsidering the problem in this paper is to 

develop a technique for thick plates which may be used in a straight- 
* 
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forward manner for solving the elasticity problem of a multilayered plate 
containing a circular hole. The practical importance of the latter 
problem lies in the fact that under unidirectional compressive loading of 
the layered medium the transverse nonnal stress on the interface is 
tensile around the hole where there is a stress concentration in a., and, 
in fact, has a (weak) power singularity at the hole boundary which greatly 
enhances the possibility of delamination failure. In this study the thick 
plate problem is formulated in a systematic manner by using the z-component 
of the Galerkin vector and that of a harmonic vector function introduced 
by Muki [2]. By using a complex eigenfunction expansion for the tractions 
on the hole boundary this formulation also leads to an infinite system of 
algebraic equations. The z-component of Muki's vector function turns out 
to be identical to the displacement potential A used by Alblas, the z- 
component of the Galerkin vector is biharmonic, and it appears that the 
potentials A, B^, and B^ of [1] constitute a special case of the functions 
ifj and Z used in this paper. In fact the characteristic equation giving the 
eigenvalues obtained in this paper is identical to that found in [1]. 

In the numerical results presented in [1] seven eigenvalues were 
computed and in the series expansion at most seven terms were considered. 
Thus, a secondary reason for reconsidering the problem is to verify and to 
complement the results given in [1] by using greater number of terms in the 
series and by providing additional results. For example, in [1] the results 
are given for v=0.25 only whereas in the perturbation problem the Poisson's 
ratio turns out to have a very significant effect on the nature as well as 
on the magnitude of the stresses. 

The problem of thick plate containing a circular hole has recently 
been studied also by Aguf and Vasil ev [3] under general axisymmetric 
boundary conditions on the plate surfaces. In [3] even though the general 
problem was also reduced to an infinite system of algebraic equations, the 
results given in the paper for a plate loaded symmetrically with respect to 
midplane of the plate were obtained essentially by a two-step successive 
approximations, each step involving the solution of a somewhat standard 
axisymmetric contact problem. In literature there are also a number of 
approximate solutions of the problem reviews of which may be found in [1-3] 
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(see, for example [4]). 


2. Formulation and Solution 

Consider an infinite plate of thickness 2h subjected to unidirectional 
stress a = a at x = ± ® (Figure 1). In cylindrical coordinates the 
stress state in the plate without a hole may be expressed as 
a a 

= -^ (1 + cos 2e) , ^ "T ’ *^zz ^ ’ 


re 


= - sin 20 


CT =0 

rz 


/^ez ^ 


(la-f) 


To solve the problem of a plate having a traction- free circular hole of 

radius a the perturbation solution of the plate obtained by using the . 

following boundary conditions must be added to that given by (1) 

a ff 

Or/a,0,z) = - (1 + cos 20) , a^g(a,0,z) = sin 20 

(2a-d) 

Orz(a,6,z) = 0 ; (a^^, 0, as r ->■ » 

The superposition of e-independent solutions given by (1) and that obtained 
from (2) gives the axisymmetric solution of a thick plate containing a 
traction-free circular hole and subjected to a uniform radial stress 
Orr = <^q/ 2 at r = «. This solution is known to be independent of z and is 
given by [1] 

a _ % /I a^> a _ % /I , a^v a _ „ 

^rr - T ■ 72"^ ’ ‘^ee ' T » '^re ‘ ° ’ 

(3) 

0^2 " ^ > '^ez ” ^ ’ *^zz " ® » Ol0f2TT) . 


Thus, the problem under consideration is reduced to a perturbation problem 
having the following boundary conditions: 


o"^(a,0,z) = - ^ cos 20 , oJ!g(a,e,z) = sin 


20 


a° 2 (a, 0 , z) = 0 


/ b b b N 

(arr.0rQ>cfr2) 


(4a-d) 


0 , as r “ . 

The complete solution of the problem of a uni directionally loaded thick 
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plate with a traction-free circular hole may then be expressed as 


a^-j{r,0,z) = ajj (r) + a^j(r,0,z) + a5^(0) (i.j = r,0,z) 

where, referring to (1) 

a o 


c _ c _ C n 

a = a = = 0 

zz rz 0z 


(5) 


( 6 ) 


The geometry of the medium and the boundary conditions suggest that 
the stress state a^^(r,0,z) , . (i,j = r,0, z) for the perturbation 

problem will have the following form; 

= Tp^(r,z) cos 20 , Ogg = TQQ(r,z) cos 28 

0^0 = 'r^e(»".z) sin 2e , cos 29 , (7a-f) 

°Qz ^ ’ °zz 

Let Z(r,0,z) and ili(r,e,z) respectively be the z-component of the Galerkin 
vector and that of the harmonic vector function introduced by Muki [2], 
The functions Z and if> satisfy the following differential equations: 


V^V^Z = 0 


(8) 




(9) 


The form of the stress state as given by (7) .and the relationship between 
the stresses and the functions Z and ^ imply that these functions may be 
expressed as [2] 


Z(r,9,z) = Zj(r,z) cos 29 
ip(r,0,z) = iJ;j(r,z) sin 20 

From (8-11) it follows that 


( 10 ) 

(11) 
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V2V|Z^ = 0 


= 0 . = 




r3r 


■ 


3 ^. 

3p" 


(12) 

(13) 


The functions Zj and defined by (10) and (11) satisfy the following 
symmetry relations: 


Zi(r,z) = - Zi(r,-z) , ^^{r,z) = 4 »j(r,-z) 

Let us now assume that.Zi and i|;^ are given by 

^^(r,z) = I f|^(r) cos B,^z , 

Z^(r,z) = I g^(r) sin a^z + z h^(r) z cos a^^z. 


(14a,b) 

(15) 

(16) 


Substituting 

from (15) and (16) into (13) and (12) we obtain 


f,(r) . 



(17) 

9„(r) = 

B„Kj(a„r) + 






(18) 

+ B„I„(a_r) + 1 
n 2 ' n ' 



h„(r) = 

C„K,(a„r) * 


(19) 

From (15-19) 

and the regularity conditions (4d) at r = » it follows 

that 

* 

B* = 0. d; = 

■ 0> C|^ ”0, (k ~ 1,2, •••) 

(20) 


Also, it can be shown that the boundary condition a^^(r,e,±h) = 0 or 
Trz(r,±h) = 0 may be expressed in the following form: 


Fi(B^,.C^, D^) F 2 (r) + D^FgCr) = 0 , (a<r«») , (21) 

where Fi, F 2 , and F^are known functions. Thus for (21) to be valid for all 
r in (a,®) one must have 
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FiCBr, C^. D^) = 0 . = 0 . (22a.b) 

By considering (20) and (22b) the functions i|»i and Zi may be expressed 
as 

'l»i(r,z) = E A,^ K2 (B|^r) cos 3,^z , (23) 


Zi(r,z) = Z^[B^ K2 (a^r) sin a^z + K2 {oL^r)z cos a^z]. (24) 

The three sets of constants A^, B^, and C^, (n = 1,2, .. .) and the 
"eigenvalues" and 3^, (n = 1,2, ...) may be obtained from the boundary 
conditions at r = a. and z ± h. Referring to (7) the conditions at 
z = ± h may be expressed as 


, ± h) = 0- , Tg2(r, ± h) = 0,(a<r<») (25) 

In terms of the functions Z and ij) the stresses in the plate are given 


± h) = 0 


by 


«b - 3 / n 2 3^ ^ 7 a. 9 3 2 X , 

'^rr ~ ~ 3r^^ 30 ^r 3r ~ r^ 


1 


-1 


32 


“re r 303z ^r ’ 3r^^ " sr^"*" 

„b _ 3 / -9 _3_ 1 _3.2 \y _3_ (Z 

00 " ^ ^ " r 3r W 30^^^ " 30 V Dr r 

<z - A 




(26a-f) 


b _ 3 


'^rz " 3 r [0-v)v2 _ ^]Z + ^ 3032 » 

-0Z=fWtO-v)V^-^]Z-^ 


303Z 

Now by substituting from (7), (10), (11), (23) and (24) into (26d-f) we 
obtain 

Tzz(»".z) = ? «n l<2(a„r){(a„ cos a„z)B„ - [(l-2v) cos a„z 


1 


n ' n 
+ »„z sin a„z]C„) 


“n '^l<“n''>‘(“n VK * sin a„z 
+ »„s cos o„z)C„) - i J^A|^6|,K2(S;,r) sin 6„z . 


(27a-c) 
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+ a„z cos o„z)C„] + A|j e,^ K^ (6|^r) s1n s,;2. 

From (25) and (27) it then follows that 

sin Bj^h = 0 , B,^ = irk/h , (28) 


(a_ COS a_h)B_ - [(l-2v) COS ah + a_h 
n n n ' n n 

(a„ sin a„h)B„ + [2v sin a„h + a„h cos 

For a non zero solution the determinant of the coefficients in (29) must 
vanish giving the following characteristic equation for a^: 

2X^ + sin 2 \^ = 0 , \ ^ » (n = 1,2, ...) . (30) 

From (30) it is seen that if is a root -X^, X^, and -X^ are also roots. 

Therefore, it is sufficient to determine the roots of (30) in the first 
quadrant of the X plane only. Also, from (29) we have 


sin = 0 , 


a h]C = 0 
n ■' n 


(29a,b) 


\ ^ + tan x„)hC„ . 


(31) 


Let us now define the following dimensionless quantities: 


\ = A|^/(aho^) 


’ C; = C^/(ahoo) , 


s = z/h 


P = r/a 


t = a/h 


(32) 


The stress coefficients x,. defined by (7) and generated by the functions 
Zi and ij)i may then be expressed as follows: 


(p»s) 


zz 


= Re E C ' [xHK,(x„pt)(tan x„ cos X„s-s sin x„s)] 
n‘‘n 2'n ' n n n 
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= - Ke J, C' [X„tK,(A„pt) . f K,(X„pt)]* 

*[X^ sin x^s + x2 tan X^ sin x^s + xjs cos x^s] 

- - ? a/ kTT Ko(kTTpt) sin ktrs , 
p k=l 

’^Q_(pjS) 2 00 I / < r • li 

= - F ^n K2( V^Jt^n V ^'n 

+ X^S COS X„s]} - Z A. kTr[kTrtKi(kTTpt) 
n n K 

+ — Ko(kTrpt)] sin kirs , 
p 

C; (-2VX2 t K 2 (X„pt) cos X„s 
-|x„ s1n X„s[X„Kj(X„pt) +^Mx„pt)] 

• ^ K 2 (x^pt) X^s Sin X^s + [2-2v + X^ tan X^]* 
*[icos X„S (X„K.(X„pt) K,(X„pt)) + 

+ K 2 (x^pt) cos 

+ Z Al cos kTTS [ -^r K 2 (kirpt) + ^ kTT Ki(kTrpt)] 

k=i k P t P 

T (p,s) 

t [x2t K„{X„pt) + ^ Kl(x„pt) + ^ K 2 (X„pt)]* 

*[x^ s sin x^s - cos x^s (2-2v + x^^ tan A^)]} 


tan x^) 
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E A. cos kirsC-^ K2(k7rpt) + — kir Ki(kiTpt)] 
(="1 ^ p 


T0(p»s) 2 


<Jo 


= - - [R^ E c' K2(X ot) + X„Ki(X Dt))* 

p e n pt n ' n ^ n " 


'[cos x^s (2-2v + tan x^) - x^ s sin x^s]}] 

- E aI cos kirs [2k2ir2 t Ko (kirpt) + Ki (ktrpt) 
k=1 , . P 

+ K2(kirpt) - k^ir^t K2(kirpt)3 

p t 


(33a-f) 


From tha characteristic equations it is seen that the problem has zero 
eigenvalues in go = 0 and Xo = 0. Therefore, for completeness some particular 
solutions must be added to li and to account for the zero eigenvalues. 

Let these solutions be of the following form: 


i|)^(r,z) = fjr) + g„(r)hjz) + mjz) 
Zi(r,z) = f*(r) + g*(r)h*(z) + m*(z). . 


(34) 

(35) 


Solutions of the form (34) and (35) satisfying the differential equations 
(12) and (13) may be expressed as 


fo(r) 

= a°, 

r2 + a“/» 

.2 

> 


9o(r) 

= 4 

fr^ . 







2b« 

ho(z) 


+ d|z + 

T^z2 

C2 

3EJ ^ 

mo(z) 


+ b°z 
2 

5 


f*(r) 

= El 

r^ 

Ear^ + 

EHr** , 


(36a-d) 
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g*(r) = G^/r^ , 

2F_ 2F, 

h*(z) = Hi + H2Z + H3Z2 + H4Z3 + 3G^ z"" + ^ » 

m*(z) = Fi + F2Z + F3Z2 + F4 z 3 . (37a-d) 

It may now be shown that imposing the boundary conditions (25), the regularity 
conditions at r = ~, and the symmetry conditions with respect to the z = 0 
plane, and requiring that the resulting stress state be non zero, (36) and 
(37) substituted into (34) and (35) give 

t|;“(r,z) = -3E(|^ + i) , . (38) 

Z°(r,z) = E(2-v) |^ + ^ + ^(l-v)z , (39) 

where E and F are arbitrary constants. Defining the dimensionless constants 



(40) 


the stress state generated by the functions and Z° may be expressed as 


T®^(p,s) 


6 F . cp' /1+y. 4 . 3vsZ 
-IT + 6E (— + 




_ 6 f‘ ISe'vsZ 




(41a-f) 


re 


(p>s) 


6F 

7 ^ 




= x° = 0 


rz zz 0z 


We now note that the stress expressions given by (33) correspond to 
nonzero eigenvalues 0^ and x^, (n=l,2, ...) and the complete solution 
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for the perturbation problem is obtained by adding (41) to (.33). Thus, for 
the perturbation problem the boundary conditions at r=a given by (4a-c) 
may be expressed as 


T^^(1»S) + T^p(l »s) = - ^ , 

TreO.s) +t;e(l*s) , (42a-c) 


T^^n.s) +Tyi,s) = 0 . 


Substituting now from (33) and (41) into (42) and expanding both sides of 
(42a) and (42b) into Fourier series in cos ttRs and (42c) into sin uks, 
we obtain the following infinite system of algebraic equations for the 


unknown constants A^, c|^, e' and f' 



Re 

00 

z 

n=l 

‘^jn 

1 

(-i)J"' 

72v 


Re 

00 

I 

n=l 

^jn 

t 


72v 

1 

a. A. + 

Re 

oo 

£ 


1 

c„ = 

0 , 

(j=l,2. 

0 J 


n=l 

jn 

n 


_ 00 
Re z 

n=l 

dnr 

or 


- 12F‘ + 

12E'(l-hj 

4. 2v\ 


Re z f^j^ C^‘ - 12 f' + 6e'(1+v + |t) = 1 • 


where 


(43) 

(44) 

(45) 

(46) 

(47) 


= - 2TrjK2(irjt) , (j=l,2, ...) , 

b-!n = - + 2K2(x„t)]Y,.„(x„ + x^ tan X„ + x^s ) , 

jn n ^ '^ n jn n n n njn 

(j»f^ ” i »2,3, .••) , 

. sin X sin x„ 

^jn “ ’ (j,n = 1,2,3, ...) 
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3 


*^jn ' ' ^ - jir (X^ - jn)^ + jiT (X^ +jii)^ 

(j.n = 1,2, 


Cj = "4TrjKi(Trjt) - ^ K2(irjt) , j-1 ,2, ...) 

dj^ = ^jn^-2vX2tK2(A„t) - {2-2v + X^ tan X^)[f (x^t) 

+ 3X^Ki(x^t) + x2tKo(x^t)]} 


+ 3X^Ki(x^t) + x2tKo(Xj^t)] , 


{n=l,2, ...; j=0,l,2, ...) 


^jn = -2t-dn^2-2v + X„ tan X^)[| K^Cx^t) + X^Ki(x^t)] 


n , ( n~l ,2, •••) , 


*^jn‘ 
sin X. 


(n=l,2, ..; j=0,l 


^on ~ ^ X 


n 


'jn = ( 


sin x_ sin x 


Xn-JTT ^n'*'^^ 


r) . (j,n = 1,2, ...) , 


cos X sin X 

<J>on = - 2 -Y-^+ 2 , (n=l,2, ...) , 

n n 

. cos X sin X cos x sin x 
‘^’jn ■ ^ “ X^-ji: (x^-jir)^ " X^+jiT (Xj^+jTT)2 ^ * 

( j »n~l , 2 , . . . 

e. = -2j2iT2tKo(Trjt) - 6jirKi(Trjt) 

0 

+ K2(7rjt)(Tr2j2t . 1^) , (j=7,2,...) 


] 

...) 


,2, ..) , 


) 

• (48) 
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The infinite system (43-47) may be solved by the "method of reduction", 
that is by truncating the series at the Nth term and by considering only 
the first N unknowns in A'.' and C.'. Note that the coefficients A.', E', and 

I 1 J J J 

F are real and are complex, (j=l,2, ...)• Thus, truncated at the Nth 
term equations (43-47) provide 3N + 2 equations to determine Ai', . . A.1 , 

I I I I N 

E, F, and the real and imaginary parts of Ci, Cj^. After determining 
these constants the stress state in the plate containing a stress-free 
circular hole and subjected to unidirectional tension at x = ±» may be 
obtained from (5). Note that the stress coefficients t (i,j=r,0,z) 

I J 

which appear in (7) are the sum of and given by (33) and (41). 

The solution and the results given in this paper are for a particular 
loading, i.e,, a = a at x = ±«. However, it is clear that the method 
is applicable to any loading condition for which the problem can be 
reduced to a perturbation problem with the tractions on the hole surface 
being the only external loads and for which these tractions can be 
expressed in terms of Fourier series in e. 

3. Results and Discussion 

As indicated before, the solution of the problem under consideration 
is given in [1]. The main reason for reconsidering the problem is there- 
fore to develop a technique which can be used in the formulation of multi- 
layered plate problems by following a somewhat more systematic approach. 

A secondary purpose of the present study is by taking advantage of the 
relative improvement in the computational capabilities since the publi- 
cation of [1], to verify and complement the results given in [1]. Thus, 
instead of seven roots of the characteristic equation in the first 
quadrant considered in [1] twenty roots have been calculated and corre- 
sponding higher number of terms have been taken into account in solving 
the related infinite system of algebraic equations (43-47). Even though 
it is very difficult to prove the regularity of the infinite system 
(43-47) theoretically, for the relative dimensions considered the con- 
vergence seems to be extremely good. Even the results given in [1] which 
are based on only seven eigenvalues do not significantly differ from those 
calculated from twenty eigenvalues. The results given in Tables 1-9 
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complement those presented in [1] in the sense that they contain the 
investigation of the influence of the Poisson's ratio which appears to 
be quite significant, the variation of the stress components other than 
and a , and the variation of the stress components in radial direction 
To show the trends of the distribution of various stresses some results 
are also given in Figures 2-7. It is important to note that in, order to 
conform to the presentation given in [1] the stresses presented in the 
tables and^in the figures are defined by (see (5)) 




(i,j = r,e,z) 


(49) 


To obtain the complete solution these stresses must be added to the 
axisymmetric stress state given by (3). 

From (43-47) it may be shown that for v=0 


A' = 0, c;. = 0, (j=l,2, ...) , E' = -l,F'=-i , (50) 

which, substituted through (41), (33) and (7) into (5), reduces to the 
following plane stress solution: 


j. "'’/I j. 3 a** 4 a^^ 

Vr " ^ ^ ~F* 


3 2 1a4 

^ (1 + -pr) - -2" (1 "pr) cos 20 


06 


2” ' ^ 


(51a-f) 


From a practical viewpoint perhaps the most important results are 

those given in Table 1 and partially displayed in Figures 2 and 3, 

namely the variation of the hoop stress on the surface of the circular 

yy 

hole. From (51) it may be observed that the "stress concentration factor" 
(ooq)»,^«/o for the plane problem has a value of 3 and is independent 
of the Poisson's ratio and the coordinate z. On the other hand Table 1 and 
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Figures 2 and 3 show that the stress concentration factor may vary with 
the location z/h, the Poisson's ratio v, and the radius-to-thickness 
ratio a/h quite considerably. The stress concentration factor is maximum 
for z=0 and is greater than the corresponding value for the plane problem. 
The table and Figure 3 show that at z=0 increases with increasing 
values of v, and approaches the plane stress and plane strain values 
{i.,e., = 2ao) as a/h approaches infinity and zero, respectively, 

becoming a maximum for some value of a/h between 0.5 and 1. For v=0.5 
this maximum may be as much as 15% greater than the corresponding plane 
stress value. 

Table 2 and Figures 4 and 5 show the calculated results for 
(-Ozz/oo cos 2e). From Table 2 and Figure 4 it may be observed that for 
a fixed value of v (-a^^/oo cos 2e) increases with the decreasing ratio 
a/h and in limit as (a/h) ->■ 0, away from the plane boundaries z=±h, it 
approaches the expected plane strain value of 2v. Also, for (a/h) -► ® 

^zz approaches zero which is the expected plane stress result. Figure 2 
implies that 

■^a^2(r,0,h> = 0 . (52) 


That this is indeed the case may be seen by considering the third stress 
equilibrium equation 


3(t t 8a„ 3 <t o 

rz ^ 1 ez ^ zz ^ _rz _ « 

3r r 30 3z r " 


(53) 


and by observing that atz=h 


or 


(r,0,h) =0 , (Tg^ (r,0,h) =0 


3r ^rz ^ * 30 ’^Qz ^ 


(54) 

(55) 


Similarly, symmetry considerations require that 
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which may also be observed form Figure 4. From Figure 5 it may be observed 
that again as (a/h) -*- « approaches zero, the expected plane stress 
solution, and for (a/h) -> 0 (-o^^/oo cos 2e) -v 2v which is the plane 
strain solution. As v and a/h vary in the intervals [0.5,0] and (0,~) 
respectively, unlike the axial stress a,, varies between the corre- 
spending plane strain and plane stress values. 

The distribution of the shear stress on the surface of the hole 
is given in Table 3 and Figure 6. Again, it may be noted that the limiting 
cases of plane stress and plane strain solutions are recovered as (a/h) -»- » 
and (a/h) ->■ 0, a., goes through a maximum in 0.5<(a/h)<l and for (z/h) = 0.8 
and increases with increasing v. 

Tables 4-9. show the spatial distribution of the stresses for a=h and 
v=0.25. It may be seen that as r » the stresses given in these tables 
added to a?, given by (3) approach the homogeneous solution (1). The 
variation oi. and with r/a is also shown in Figure 7. 

ou ZZ 

The first twenty roots of the characteristic equation (30) are 
given in Table 10. They seem to agree with the first seven roots 
calculated by Alblas in nearly all significant digits. 
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TABLE 2 


Distribution of the normalized axial stress boundary r=a 





























TABLE 3 


Distribution of the normalized shear stress -o^^/ao sin 2e at the hole boundary r=a 
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Variation of cos 20 with r/a, =1 , v = 0.25) 
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Figure 2. Thickness distribution of the normalized hoop stress 

-a^g/ao cos 20 at hole boundary r = a for v = 0.25 and for 
various values of a/h. 



□ / h 


Figure 3. Variation of the normalized hoop stress wiht a/h calculated 
at z = 0, r = a. 




Figure 4. Thickness distribution of the normalized axial stress 

cos 29 at the hole boundary r = a for v =0.25 
and various values of a/h. 





Figure 6. Variation of the normalized shear stress -a. /oo sin 2e 

QZ 

with a/h calculated at r = a and z = h/2. 



Radial distribution of normalized hoop and axial stresses 
calculated at z = 0 plane for v = 0.25 and a = h. 


Figure 7. 
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